CoUisional relaxation of Feshbach molecules and three-body recombination 
in ^'Rb Bose-Einstein condensates 
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We predict the resonance enhanced magnetic field dependence of atom-dimer relaxation and three-body re- 
combination rates in a '^'Rb Bose-Einstein condensate (BEC) close to 1007 G. Our exact treatments of three- 
particle scattering explicitly include the dependence of the interactions on the atomic Zeeman levels. The 
Feshbach resonance distorts the entire diatomic energy spectrum causing interferences in both loss phenomena. 
Our two independent experiments confirm the predicted recombination loss over a range of rate constants that 
spans four orders of magnitude. 
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Few-body collisions determine the lifetimes of cold gases 
of atoms and diatomic molecules. Recent experiments on 
identical Bose atoms Uj, |31 suggest that magnetically tun- 
able interactions Ij, 01 in combination with such scattering 
phenomena could allow association of atomic clusters at rest 
fSD and confirmation of long standing predictions of quantum 
physics II2]. Destructive interferences of inelastic three-body 
collisions have paved the way for the Bose-Einstein conden- 
sation of '■'''Cs [2]. Scattering processes involving three atoms 
are known as three-body recombination and atom-dimer col- 
lisions. Three-body recombination refers to the threshold-less 
transition from initially unbound atoms to a dimer molecule 
and a remnant atom in accordance with energy and momen- 
tum conservation. All three atoms are lost from a cold gas. 
Previous theoretical studies relied upon model calculations 
131 or analytic estimates |8], suggesting general trends for 
the behaviour of recombination rates in the limit of large di- 
atomic scattering lengths. Several ab initio approaches to cold 
atom-dimer collisions |9] are presently limited to compara- 
tively tightly bound diatomic states with energies up to about 
-A:b X 1 K (fefi = 1-38 X lO^^^ J/K). The energies of dimers 
produced via magnetically tunable interactions are, however, 
typically on the order of -k^ x 1 mK, right below the dissoci- 
ation threshold. 

In this paper we predict, for the first time, the atom-dimer 
collision rates involving such Feshbach molecules, consist- 
ing here of **^Rb atoms, without any adjustable parameters. 
Our exact approach to the three-body Schrodinger equation 
explicitly accounts for the coupling between the atomic Zee- 
man states. This allows us to properly describe the magnetic 
field dependent distortion of the diatomic bound state energies 
I IC, Uin caused by the weakly coupled Feshbach resonance 
1 12l[l3ll . We predict constructive and destructive interferences 
in both atom-dimer relaxation and three-body recombination 
rates which we interpret in terms of resonance phenomena as- 
sociated with meta-stable trimer or bound dimer levels. Our 
approach has been tested by two independent experiments on 
three-body recombination in a '^^Rb BEC over a range of rate 
constants that spans four orders of magnitude. 

The three-body Hamiltonian, //jb, is comprised of the ki- 



netic energies of the atoms and their pairwise interactions, 
i.e. //3B - Hct + V\ + V2 + V-},. Here V\ is the potential asso- 
ciated with the interaction of the atom pair (2, 3), whereas V2 
and Vt, follow from cyclic permutations of the indices. We ne- 
glect short range genuinely three-body interactions as they are 
not expected to significantly affect resonance enhanced colli- 
sions. In the barycentric frame the total kinetic energy, //q, 
can be divided into contributions from the relative motion of 
an atom pair and the motion of the third atom with respect to 
the centre of mass of the pair. Given the associated spatial 
Jacobi coordinates, r and p, respectively, such a separation 
yields: 



^o = |bg)|--V^-V?|(bg|. 



(1) 



Here m is the atomic mass and |bg) indicates a product of (/ = 
1 , mj - \) Zeeman states 1 11 ] in which the three atoms of the 
"^^Rb BEC are prepared. Depending on the context, we refer to 
either a product of three or two of such atomic ground states 
as the entrance-channel spin configuration and choose the zero 
of energy at its threshold for dissociation into unbound atoms. 
To describe their initial and final configurations, we label 
each one of the colliding atoms by a Greek index a - 1, 2, 3. 
The pairwise interactions, y„, determine the properties of cold 
diatomic collisions and highly excited molecular bound states, 
as well as their dependence on the magnetic field strength, B. 
In an idealised treatment, the observed resonant enhancement 
of the binary scattering length in the vicinity of B() = 1007.4G 
Ill2i [l3l1 arises from the near degeneracy of a single meta- 
stable diatomic energy level with the entrance-channel disso- 
ciation threshold. This bare Feshbach resonance level, |0res), 
refers to a closed scattering channel comprised of Zeeman 
states from the excited / = 2 le vel jllll . Accordingly, the 
interactions can be represented by fill] 



V«=|bg)y,^^(bg| + W„|0,-es).„( 



\<t>Ki)aa{(l>K&\Wa + \(pKs)a E^^iB) a{ 



(2) 



/bg 



Here Vj^ and |(^res)ff are the diatomic entrance-channel poten- 
tial and resonance level, respectively, whereas the third atom 



with index a is in the Zeeman ground state and plays the role 
of a spectator. Similarly, ]¥„ describes the pairwise inter- 
channel spin exchange coupling. The resonance state energy, 
Eies, depends on B through Er^siB) - fire^iB - Bres)- Here 
//res - h X 4.2 MHz/G is the difference in rnagnetic moments 
of closed- and entrance-channel atom pairs 1 1011 . and Bies indi- 
cates the point of degeneracy of Ziies with the entrance-channel 
dissociation threshold. As the resonance state is weakly cou- 
pled, the shift Bo - Bres is negligible lUlll . Due to the absence 
of overlap between Zeeman states constituting the diatomic 
entrance (/ = 1) and closed (/ = 2) channels, all resonance 
states associated with different arrangements of the atoms are 

orthogonal: „(0i-esl<;*res>/J = Sap- 

Our exact solutions of the three-body Schrodinger equation 
are based on the Alt, Grassberger, and Sandhas (AGS) tech- 
nique [14]. Accordingly, we introduce the complete Green's 
function, Gisiz) = (z - /^sb) ', which characterises all bound 
and continuum levels of Hj,^. Here z = Ei + iO is a com- 
plex variable indicating that the physical continuum energy 
of incoming particles, Ei, is approached from the upper half 
of the complex plane. Given the arrangement-channel and 
free Green's functions, Ga{z) - {z - Hq - VaY^ and Go(z) - 
{z - Hoy^, respectively, the AGS transition matrices, Uapiz), 
are defined by implicit relations (a,/3 =1,2, 3): 



10" 



G3b(z) = dapGiiiz) + Ga(z)Ua/l(z)G,3(z). 



(3) 



In accordance with Ref. |14], these transition matrices de- 
termine, e.g., the Bose-symmetric scattering amplitudes for 
atom-dimer elastic or relaxation collisions by 



/(qf,qi) 



in mh 



Y^a{({i,<li\Uap{zMJb)p- (4) 



»,/3=l 



Here |qi,0j,)/? indicates an incoming dimer in the coupled- 
channels state |0j^) of energy E'^, whereas the third atom with 
index /? has the momentum qi with respect to the molecu- 
lar centre of mass, and is described by (plqi) = exp(/qi ■ 
plK)l{2nKfl^. Consequently, the initial state fulfils {Hq + 
V/j)lq„<^i,)/j = £,lq„<^i,)/j with E, = 2>qll{Am) + E[. Simi- 
larly, iQf , (^{^)„ describes the scattered atom and dimer, whose 
final energy, Ef - 3q^/{4m) + E^, equals Ei. The transition 
matrices in Eq. (|4]l are determined by the AGS equations il4ll : 

3 
Ua0(z) = (1 - 5„^)Go'(z) + Yj VyGy{z)Uy,3(z). (5) 



Based on Eq. Q, the upper panel of Fig. [T] illustrates the 
rate constants, /Tad = o-(qf)qf/{2m/3), for elastic scatter- 
ing {E*^ - E'^ - £■_!) as well as relaxation (is^ < E'^ - 
E-i) of **^Rb2 Feshbach molecules. Here o-{qf) denotes the 
(in)elastic i-wave cross section associated with the amplitude 
of Eq. dUl. Our implementation of Eq. Q follows the ap- 
proach of Ref. f 15] . Consequently, the interactions properly 
describe: the measured B-dependences of the diatomic scat- 
tering length [13, 16], the binding energies of the two highest 
excited vibrational '^^Rb2 states 1.10,] . their long-range wave 
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FIG. 1: (Colour online) Upper panel: Elastic (dashed curve) and 
inelastic (solid curve) atom-dimer scattering rates, K.^^, associated 
with a collision energy of 3q-/(4m) = k^ x 10 nK (k^ x 1 mK= 
h X 20.8 MHz) versus B on the low-field side of Bq. Lower panel: 
Energies of the Feshbach molecular (dashed curve) and second high- 
est excited vibrational (dotted curve) dimer states, £_i and £'_2, re- 
spectively. The dot-dashed and solid curves indicate the energies of 
isotropic meta-stable **''Rb3 states. The highest excited trimer level 
emerges at about 1003.7 G (arrow) in a zero energy resonance (in- 
set) for collisions between atoms and dimers in the second highest 
excited vibrational state. 



functions as well as their entrance- and closed-channel spin 
admixtures |11]. The next vibrational dimer level is far de- 
tuned from the dissociation threshold by about 600 MHz 1 15[] 
and will be neglected. Our predicted collisional relaxation 
rate constants are typically on the order of 10""' cm^/s. Such 
magnitudes have been confirmed by experiments on the sta- 
bility of ^^Rb2, yielding /Tad = 2 x lO^'^cm^s at 1005.8 G 
[17] in good agreement with Fig.[Tl as well as for ^^Na2 Fesh- 
bach molecules [18]. Both species are associated with weakly 
coupled, closed-channel dominated [ 1 1[, resonances. The fast 
atom-dimer decay tends to exceed elastic scattering, except 
for a region centred at about 1001 G, close to the avoided 
crossing of E-i and E-2 at 1001.7 G illOtl in the lower panel 
of Fig. \T\ We expect relaxation into more tightly bound dimer 
states to partly fill in the gap in the inelastic rate constants. 

Three-body recombination in a BEC is well described in 
terms of the transition from a continuum momentum state, 
|0,bg), associated with separated atoms at rest, to a dimer 
and a remnant atom. Whereas the initial state thus fulfils 
//o|0,bg) - 0, the final states are of the same nature as those 



in atom-dimer relaxation and elastic collisions. In accordance 
with Ref. 1,14.1 . the associated transition matrix is given by 
Ua,0=o(z) of Eq. (O with z = /O. For a definite target dimer 
state 10^) of energy is^ this matrix gives the loss rate constant 
to be: 

^3 . W2^^)^ j^ r^Q|„(q,4|f/.o(/0)|0,bg)f . (6) 
(2=1 

Here dQ. refers to the angular component of qt whose modu- 
lus, qf, is determined by 3q^/{4m) + E^ - 0. The total loss 
rate constant, .^3°', is found by summation of Eq. ^ over all 
target states. Consequently, the number of condensate atoms, 
A^, decays in accordance with 



A^ = -K'°\n-)N/6. 



(7) 



Here (n^) is the mean square density and the factor 1/6 ac- 
counts for the coherent nature of the gas 11911 . The transition 
matrix, t/g-pCiO), can be inferred from the solutions of the AGS 
equations lfl4ll using jS = in Eq. ^. 
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FIG. 2: (Colour online) Three-body recombination loss rate con- 
stant, Ki, versus B. Solid curves refer to the total loss rate constant, 
A""', whereas dashed and dot-dashed curves are its contributions from 
the final Feshbach molecular and second highest excited vibrational 
*^Rb2 states, respectively. The inset shows an enlargement close to 
the diatomic resonance whose position, 60, is indicated by the dotted 
line. 

Figure |2] gives an overview of the predicted total loss rate 
constants, /T^"', as well as their contributions from the different 
dimer target states. Recombination into Feshbach molecules 
dominates Kf^ typically by an order of magnitude. Similarly 
to the atom-dimer collision rate constants of Fig. [1] the loss 
rate constant /^°' shows interference minima which may be 
attributed to B-dependent changes in the bound states of //3B. 
Such isotropic trimer levels are indicated by the solid and dot- 
dashed curves in the lower panel of Fig. [T] and were deter- 
mined using the Faddeev approach of Ref. |5]. The minimum 
of /T™ at about 1004G in Fig.|2]occurs close to an atom-dimer 
zero energy resonance indicated by the arrow in Fig. [1] At 
the resonant magnetic field strength the trimer level becomes 



degenerate with the second highest excited diatomic vibra- 
tional state where it can decay into a dimer and an atom in 
accordance with energy conservation. Both trimer levels of 
Fig. [T] are unstable with respect to such decay processes in- 
volving more tightly bound dimer states with energies below 
-/ix 600 MHz. In the inset of Fig.|2]the interference minimum 
below Bo and the maximum just above Bq can be attributed 
respectively to atom-dimer and three-body zero energy res- 
onances associated with meta-stable Thomas-Efimov trimer 
levels 12. SH). Similarly, the oscillations in the atom-dimer 
collision rates close to Bq in Fig.[T]result from Efimov's effect 
[6], but their narrow magnetic field range makes this physics 
largely inaccessible to experiments on ^^Rb. The minima of 
Kf' at about 1004 G and 1008.5 G in Fig.|2l however, could 
be resolved. 
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FIG. 3: (Colour online) Comparison between measured (circles and 
squares) and predicted (solid curve) total three-body recombination 
loss rate constants, Kf\ Error bars refer to statistical uncertainties 
only. The inset shows the binary scattering length, a(B) = abg[l - 
AB/(B - Bo)], on the low-field (solid curve) and high-field (dashed 
curve) sides of the diatomic zero energy resonance (vertical dotted 
line). The horizontal dotted line indicates the background scattering 
length (13, [I6] flbg = 100.5 flBohi (flBoin = 0.0529 nm) whereas the 
resonance width is AS = 0.21 G. All measurements refer lo B < Bq 
where the BEC is stable with respect to collapse. 

Figure [3] compares the predicted Kf'^ as a function of B 
with two independent experiments. In accordance with the 
inset, the measurements cover a range of diatomic scatter- 
ing lengths from about 100 to lOOOaBohr within which Kf^ 
varies over four orders of magnitude. In the Oxford exper- 
iment, indicated by circles in Fig. [3] **^Rb atoms are loaded 
from a BEC into about 60 wells of a one dimensional optical 
lattice. The lattice is formed from a retro-reflected, vertical, 
850 nm laser beam giving rise to individual wells with trap- 
ping frequencies (47,47, 1.75 x 10^) Hz. The atoms are then 
prepared in the (/ = 1 , nif = 1 ) Zeeman ground state and a 
homogeneous magnetic field is ramped to a value B just be- 
low Bo- At this stage the atom number is A^o = 3 x 10^ and 
the mean square density is («q) - 2x 10^^ cm"^. This is ob- 
tained by calculating the mean square density in each well 
using the BEC wave function from Ref. [20] and averaging 



over the populated sites. After a hold time of Af = 100 ms, 
the magnetic field is ramped off and the final atom number 
A^f is measured. Assuming that coherence is maintained and 
that the volume of the cloud used to calculate the density n 
in Eq. ^ remains unchanged, the loss rate constant is ex- 
tracted from ^3°' = 6 [(No/Nff - l] /(2<n^)Af). Uncertainties 
in A^o prevent measurements of Kf^ below ~ 5 x 10"^^ cm^s. 
We have confirmed this analysis by dynamical simulations of 
the experimental sequence using the non-linear Schrodinger 
equation with a three-body loss term. 

In the Garching experiment, indicated by squares in Fig. [3] 
a BEC of 5 X 10^ atoms is loaded into a crossed-beam opti- 
cal dipole trap |16]. After preparing the gas far away from 
Bo the magnetic field is ramped and held at its final value, B. 
For B < 1007.2 G the size of the cloud changes adiabatically 
during the experimental sequence. Using the Thomas-Fermi 
prediction for the mean square density, /Tj"' is obtained from 
a fit of Eq. © to the loss data |21]. For B > 1007.2 G the de- 
pletion is too fast for the BEC to adiabatically adapt its size. 
Consequently, we jump the magnetic field from the initial to 
its final value and measure the loss only on time scales suffi- 
ciently short for the cloud to keep its initial size. Systematic 
errors arise mainly from the calibration of the trap frequencies 
and the atom number, leading to an estimated uncertainty of a 



factor of 3 in /Tl". Statistical errors are negligible in compari- 
son. The accuracy of magnetic field calibration is 30 mG. We 
believe that our B-fields are sufficiently far from Bq for atom 
loss caused by collisional avalanches |22] to be negligible. As 
the BEC dynamics in the crossed-beam dipole trap is compar- 
atively simple, the associated overall systematic uncertainty is 
smaller than in the optical lattice setup. The agreement within 
uncertainties between predicted and measured K!^^ of both ex- 
periments in Fig. [3] confirms our approach. 

We have shown how spin-dependent potentials can be in- 
cluded in a practical and exact treatment of three-body scatter- 
ing phenomena. Destructive interferences predicted in atom- 
dimer relaxation rates allow elastic scattering of atoms from 
generally rather unstable ^^Rb2 Feshbach molecules. The pre- 
dicted magnetic field dependent suppressions of three-body 
recombination can be used to increase the lifetimes of cold 
gases [2J. While we have confirmed our approach by exper- 
iments on ^^Rb condensates, it is applicable to all species, 
bosons and fermions, subj ect to closed- or entrance-channel 
dominated resonances flltl. 
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